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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

APPLICATION OF STATISTICAL THEORY TO BEAM-RIDER 
GUIDANCE IN THE PRESENCE OF NOISE. 

II - MODIFIED WIENER FILTER THEORY 
By Elwood C. Stewart 


SUMMARY 




A study has teen made of the application of Newton’s modification 

I of the Wiener filter theory to the optimization of a beam-rider guidance 
system operating in the presence of glint noise. Target and missile 
motions are assumed to be coplanar. 


- The theory has been applied using a typical variable -incidence mis- 
sile by placing a realistic restriction on the mean-square surface deflec- 
tion so that the system operation is confined to the linear range. The 
transfer functions of the optimum guidance system are derived and an 
example missile-control system is synthesized. It is shown that the 
min i mum attainable error corresponding to a realistic control restriction 
is close to that for the Wiener theory. Performance of the system versus 
noise magnitude is given. 


It is shown that the most critical saturable quantity is control 
deflection since constraining control deflections to realistic values 
can prevent saturation of other important quantities. Servo energy 
requirements are also greatly reduced in comparison with systems in 
which saturation is allowed to occur. 





INTRODUCTION 


Noise effects in missile-guidance systems impose one of the most 
serious limitations on the effectiveness of a missile. Noise signals in 
general have the effect of adding to the system false information which 
cannot be distinguished from true information. Thus the missile responds 
to the noise signals as well as to the true signal; the miss distance is 
thereby increased. Since many sources of noise (such as angular scintil- 
lation) are inherent in the physical mode of operation and cannot be 
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removed, it is important that missile systems he designed to minimize 
the miss even though the noise is present. Statistical filter theory is 
useful in performing this minimization. 

One of the most valuable theories has been developed by Wiener 
(refs . 1 and 2 ) . A previous report has considered the application of this 
theory to the minimization of the effects of radar glint noise in a beam- 
rider guidance system (ref. 3) • It was found there that the optimum per- 
formance specified by the theory coiold not be obtained because of certain 
practical restrictions, as for example limiting of control deflection. 
Further study showed, however, that even in the presence of limiting, the 
optimum performance could be approached by the addition of network filter- 
ing to the guidance system. On the other hand, the optimum form for the 
added filtering could not be determined from this study since the system 
design was based on a theory which did not take into account the finite 
range of operation of the saturable quantities. In the present report an 
attempt is made to overcome this deficiency by determining the optimum 
transfer function under this last restriction. The theory for such an 
approach was made available in a paper by G. C. Newton (ref. 4) on a modi- 
fication of the Wiener theory. 

In the application of this theory to the missile guidance problem it 
is necessary to make certain assumptions. The main assumption made in 
reference 3, that is, the assianption that the target and missile move in 
the same plane, also will be made here. Other assumptions, such as those 
relating to the class of target maneuvers and noise, are discussed in the 


SYMBOLS 


transfer ftmction of optimum compensating network 
transfer function of the fixed network 

noise magnitude or zero frequency spectral density, ft^ /radian/sec 
time constant of the noise spectrum shaping filter, sec , 
optimum closed-loop transfer function 
acceleration of target maneuver, ft/sec^ 

twice the average switching rate of target acceleration, l/sec 
target displacement from a space reference, ft 

apparent target displacement from true target center due to noise, 
ft 
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missile displacement from a space reference, ft 
$id cross -spectral density of the input signal with the desired output 

spectral density of- the input signal 

spectral density of noise displacement yjj, ft^/radian/sec 
spectral density of target displacement y^i, ft^/radian/sec 
error between target and missile position, y^i - y^^, ft 
component of error e due to target motion, ft 
component of error e due to noise, ft 
Lagrangian multiplier 
control- surface deflection, radians 
angle of yaw, radians 
optimum open-loop transfer function 

open-loop transfer function of system approximation to [Iq 
angular frequency, radians /sec 
complex conjugate of ( ) 


GENERAL CONSIDERATIONS 


Glxnt noise is a term that is used to describe a shift in the appar- 
ent target location as determined hy a radar. It is due basically to the 
variable reflection characteristics of aircraft targets and arises from 
the relative movement of the various reflecting surfaces. Since the radar 
utilizes the reflected signal to determine target location, variations in 
the reflected signal are interpreted by the radar as shifts in the target 
location. Of the many sources of noise which may exist in a missile- 
control system, glint noise is one of the most serious since it has a 
physical origin which cannot be eliminated, imposing a fundamental limi- 
tation on missile effectiveness. The situation in regard to glint noise 
is illustrated in figure l(a) iirtiere the true target position is indicated 
as y^ and the glint noise is represented by the displacement yjg* 
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The present report is restricted to a two-dimensional study in which 
the target and missile move in a horizontal plane. ^ The guidance system 
is considered to he of the heam-rider type illustrated in figures 1(a) and 
1(h) . It should, he noted here that displacements are referred to a 
fixed space reference. The ^general problem of the guidance system is to 
make the missile position coincide as closely as possible with the true 
target position yip. The difference yip - yj^ is indicated on this figure 
by the error €, which obviously should be minimized. The criterion of 
merit which will be used here is the conventional mean-square time average 
of the error. This criterion is particularly appropriate in the case of 
the beam-rider system because the target-to-missile range is not normally 
transmitted to the missile; since the missile never knows when the target 
will be reached, it is reasonable to minimize the error for all values 
of time . 

The design of a system normally depends on the inputs to be encoun- ' 
tered, in this case the target motion and the noise. Because of their 
random nature neither of these quantities can be conveniently described 
explicitly as functions of time, and statistical descriptions are more 
suitable. In the following paragraphs the target motion and noise are 
briefly described. 

The glint noise can be defined by specifying both the power spectrum 
or power spectral density as a function of frequency, and the amplitude 
distribution. Many measurements have been made on the characteristics of 
glint noise. (For a brief bibliography see ref. 3*) Although these meas- 
urements are somewhat complicated and uncertain, it is generally found 
that the amplitude distribution is approximately Gaussian and that the 
spectral density can be adequately represented by 

" Tj^2(j2 + 1 (1) 

An examination of glint noise spectra indicates that the break point 
p/2jtTjj) is generally on the order of several cycles per second and as 
in reference 3 will be taken to be 6 cps, corresponding to Tw = O.0265 
second. The magnitude of the spectrum, W, of equation (l) depends on 
factors such as target size and target aspect so that the guidance system . 
is generally forced to operate over a wide range of magnitudes . This 
range imy extend from 7 ft^/radian/sec for small targets up to around 
30 ft^ /radian/sec for large bombers. For this range of magnitudes it 
might be thought that it would be necessary to optimize the system for 
each noise magnitude. As shown in reference 3> however, this was not the 
case for the method used therein since near-optimum results were obtained 
by optimizing only for a mid-range value called the design value. It is 

^Ibe complete three-dimensional problem would require a more complex' 
analysis than used herein. Possibly either the present theory or Wiener’s 
th^ry for multiple time series (ref, l) could be applied to this case. 
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not unreasonable to believe tbat the same will hold for the method of the 
present report. As in reference 3> ^ design value of 15 ft^/radian/sec 
will be used here. 


Consider now the target maneuver. It is difficult to say exactly 
how a target will maneuver when under attack. A reasonable situation 
might be one in which the target maneuvers in some random manner to avoid 
being hit. As in reference 3 ) it is assumed that the target maneuvers 
laterally with maximum acceleration alternately in opposite directions. 
The duration of each acceleration is a random function determined by some 
distribution. A reasonable distribution which leads to an_easily handled 
spectral density is the Poisson distribution [ (l/T.)exp(-T/T) ], where T 
represents the time and T the average time between switches of the 
acceleration. The spectral density of the target acceleration for this 
maneuver is described by 

ka^ 

^T + k^) 


Here the quantity a represents the magnitude of the target acceleration 
normal_to the beam, and k is twice the average switching rate or 
k = 2/T. The spectral density of the target displacement is then given 


■by 


^T = ~ 


ka^ 

+ k^) 


( 2 ) 


Because of the w"* in the denominator it would appear that there is a 
problem here concerning the existence of this spectral density. However, 
it can be shown that it is possible to use this representation for pur- 
poses of computations (see Appendix A). For the typical tail-chase maneu- 
ver to be used in a later example the target is assumed to maneuver with 
±1 g acceleration at an average period of 5 seconds, which gives 
a = 32*2 ft/sec^ and k = 0.4 switch/sec. 


It should be pointed out that an optimum design based on the above 
target maneuver will operate efficiently against this class of maneuver 
as a whole, although not necessarily efficiently against one particular 
target maneuver such as a single target turn. However, as indicated in 
reference 3> this design will be essentially optimum for the single-turn 
maneuver as well. 
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analysis/ RESULTS^ AND DISCUSSION 
The Modified Wiener Theory 


The primary objective of the modified Wiener theory is to arrive at 
an optimum transfer function in 'which the effects of certain saturating 
elements are considered. Ideally it would be desirable to determine this 
optimum transfer function ■without any assumptions about whether or not the 
system would operate in these saturation regions. However, no such theory 
is available at present. The approach considered herein is one in which 
the transfer. functions of certain available elements are chosen so that 
saturation does not occur. Thus the system is forced to remain linear 
and linear methods are immediately applicable. Although it- is not appar- 
ent at first glance, there are certain indications that this approach may 
lead to desirable results. First, it appears intuitively that saturation 
is undesirable because of the loss of both intelligence and ability to 
control the missile when in the saturated regions. Second, it is indi- 
cated in reference 3 that the undesirable effects of saturation can be 
reduced by the addition of filtering without undue increase in the error; 
this method in effect tends to keep the operation of the system in the 
linear range. These are only indications, however, and the real justi- 
fication in the completely linear approach depends on the resultant error 
performance which can only be established by investigation. 

The theory to be used in this approach follows Newton's modification 
of the Wiener filter theory. To describe the essentials of the theory it 
will be convenient to refer to figure 2 where the block diagram corre- 
sponding to the general filter problem is given. It should be noted that 
the notation here corresponds to that in reference In the general 
situation there are usually certain elements that are given and may not 
be altered in the design of the system. These elements are denoted by 
the transfer function . in the figure. In particular cases the fixed 
elements might represent missile aerodynamics, servo motors, and so forth, 
depending on the application. To describe the limiting or saturating 
quantities it is necessary to express them in terms of the fixed elements 
and either the input or output of the system. Which choice is made de- 
’ pends on the location of the fixed elements - if the output element is 
fixed, as is usually the case, it is desirable to express the saturating 
quantities in terms of the output. Figure 2 illustrates a situation 
wherein the output element is fixed. At this point it is necessary to 
know the saturating quantities. For this purpose we introduce a ficti- 
tious situation. Since in the, general case the saturating quantities are 
different from the input to the fixed network, the input to in 

figure 2 is imagined to be fed through certain transfer functions Hsn 
which act as recording elements and whose outputs are the saturating 
quantities 0sn* I* should be pointed out that according to this formu- 
lation of the problem the outputs of the saturating elements do not feed 
into the system but are used merely for the purpose of evaluation. From 
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tlie above discussion, then, it is clear that by specifying both the fixed 
and saturating networks, the saturating quantities can be related to the 
output. The remainder of the system, indicated in figure 2 by the compen- 
sating network is unknown and must be determined according to a 

specified criterion of merit. 


By way of comparison, the original Wiener theory considers only the 
problem of choosing the over-all transfer function so as to result in 
minimum mean-square error, 06, between the actual output, 9 q > and the 
desired output, 0^* The modified theory, however, considers the minimi- 
zation of this same mean-square error with a side restriction on the 
available range of operation of the limiting quantities; for mathematical 
reasons it is convenient to consider this restriction to be in the form 
of a mean-square limitation. According to the method of Lagrange (ref. 5) 
this means that the quantity to be minimized is of the form 

W 

^ Pn®sn (3) 

n=i 

where the. p^’s represent Lagrangian multipliers which must be chosen 
properly to restrict the operation of the saturable quantities 0gn to 
very nearly the linear range. 


In most cases it is possible to control saturation by an appropriate 
choice of the mean-square value of the saturating quantities Qgn. In the 
simple case when a quantity has a Gaussian distribution of amplitudes the 
relation between the mean-square value of this quantity and the probability 
of its limiting is well known and particularly simple. The distribution of 
the saturating quantity in the actual case may not be Gaussian for either 
of two reasons: First the input may not be Gaussian or, second, there may 

be several saturating quantities which interact. On the other hand, even 
in these cases the distribution is frequently near enough to Gaussian that 
saturation can be readily controlled by an appropriate choice of, mean- 
square value. / 


The derivation of the solution for the optiminn compensating network 
which minimizes the expression ( 3 ) is beyond the scope of this report, 
but as shown in reference an expression for this network is 


Hco(io)) = 


2«A+(u)) Jo 




Hf(ia)$j^^(a)e 
5 A"(a) 


iat 


da dt 


where 


N 

A(u) = (iw)Hf (iw) PnHgn(iw)H3n(iu) 

n=i 


> ( 1 ^) 
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In this equation the <I> quantities represent certain combinations of 
target motion and noise spectral densities (see Appendix A), refers 
to the fixed network transfer function, Hsn refers to the transfer func- 
tion of the nth saturating component, and Pj^ represents the correspond- 
ing Lagrangian multiplier. Of these quantities, all will nonnally he 
known with the exception of the Pjj. This means that only the fonn of 
Hco is known from this equation since the numerical constants depend on 
the values of The p^? however, are implicitly related (ref. !»•) to 

the mean-square restrictions on 0sn following equation 


pOO 

= / HgQ(iw)HQQ(iw)Hgj^(iu))Hgj^(i(i))$j^j^((j)d(A) (5) 

W _oo 


Here it is to he noted that the integrand in equation (5) is dependent 
on the Pi^P2> • • • virtue of the dependence of Hco these quan- 

tities (eq. (i)-)). It can he shown that in many cases of interest the pn*s 
are uniquely determined hy the equations ( 4 ) and ( 5 ) "vdien a suitable mean- 

square value, 9sn^ , is chosen. This quantity is to he chosen so as to 
limit the probability of saturation to a small value according to a cri- 
terion given later. In most practical applications the complexity of the 
operations involved in solving these equations does not permit the general 
solution to he obtained explicitly in terms of the input characteristics 
and fixed networks. In spite of this a great deal can he learned from 
certain numerical solutions, as will he shown. 


Application of Modified Wiener Theory to the 
Beam-Rider Guidance System 


The validity of the application of the modified Wiener theory to the 
beam-rider guidance problem depends on whether or not several restrictions 
inherent in the theory are met. In general, they are no stricter than are 
those of the unmodified Wiener theory. First, the input quantities, target 
motion and noise, must he stationary random series (see ref. 6 for a 
detailed definition) and defined hy corresponding power spectra which are 
continuous. It is generally believed that displacements at the target are 
approximately stationary random series (ref. 3 ) • Since the beam-rider 
system operates from these displacements, the inputs to the beam-rider 
system are also stationary random series. Second, the transfer function 
of the system must have constant coefficients. In general, .however, the 
kinematic loop of guidance systems involves a time-variable range factor. 

In particular, for the beam-rider system the time-variable factor is the 
ratio of the launcher -to -missile and the launcher-to-target ranges. On 
the other hand, because the miss is determined primarily by what happens 
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near the end of flight during ^ich the variation in this ratio is small, 
it is reasonable to assume that the requirement of constant coefficients 
is approximately met. Third, the transfer functions of the fixed and 
, saturating components must he known. 

Simplifications . - There are many quantities which may saturate in the 
heam-rider system; the most important are (l) control motion, (2) rate of 
control motion, and (3) radar receiver voltages. To apply the theory to 
all these nonlinearities simultaneously would he a tremendous task. How- 
ever, there are indications that only one limiting quantity need he con- 
sidered in the present problem. 

One indication is that of the three types of limiting not all are of 
equal severity so that one type tends to predominate. Intuitively, it 
might he felt that control-motion limiting is the most serious type and 
that placing a restriction on the control motion will satisfactorily reduce 
the other types of limiting. Of course the validity of this presimiption 
would require verification. 

Another indication is that of the three types of limiting one is more 
fundamental in that it is more difficult to remove than the others. It 
will he noted that items (2) and (3) above are both control-system limita- 
tions while (1) is essentially an aerodynamic or mechanical limitation. 
Since control systems are relatively versatile and can he designed to cover 
a wide range of characteristics, it appears that control-motion limiting 
is the most fundamental to the problem. 

From these indications it appears reasonable to consider only control- 
motion limiting. The validity of the assumption that reducing this one 
type of limiting will reduce the other types sufficiently will then he 
verified. In terms of figure 2 this simplification means that only one 
saturating-component transfer function, Hg, is involved. If the fixed 
network is chosen to he the aerodynamic transfer function so that the 
control motion is the input to the fixed network, then Hg equals unity. 
Thus, a block diagram more suited to the specific case of the beam-rider 
system would appear as in figure 3- I’t will he noted in figure 3 that the 
input quantity called 9 ± in figure 2 has been split into the two parts 
which exist in the actual case, ^target motion yiji and noise y^^. Like- 
wise the desired output 0^ of figure 2 is represented in figure 3 "by the 
target motion y,p. The box Hco^ then, is that portion of the guidance 
system which is to be chosen to minimize the expression (3). It may be 
noted that in the form given the diagram does not bear a direct resemblance 
to the form of an actual guidance system since the latter will involve 
feedback loops. The form shown, however, typifies that of the general fil- 
ter problem and is better suited to calculation. The solution of this 
problem can be readily converted to that of the actual beam-rider system. 
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A few additional remarks are appropriate concerning the two quanti- 
ties and Oil in equation (!<■) . If $ii is defined as the autocor- 

relation function of the input in figure 2, then 

1 ' 

cPii(T) = lim ^ / 0i(t)0i(t + T)dt 

T^od 

In terms of the input of figure 3, this becomes 

, pT 

cpii(T) = lim ^ / [y^Ct) + yN(t)][yT(t + t) + yjj(t + T)]dt 

It is normally assumed that there is no correlation between target motion 
and noise so that cpitn(t) and cPj^(t) are zero. It follows that the Fourier 

transform of 

= $,p(u) + %(w) 

Since in the present case the target motion and the desired output are 
identical it is easy to show in a simila,r fashion that.. 

^id(‘*>) = 

General solution with control motion restriction .- By virtue of the 
simplifications discussed in the previous paragraphs the application of 
the theory to a beam-rider-type guidance system becomes feasible . For the 
case illustrated by figure 3# optimum compensating network, ®co^ 
given by 

Hco(iw) = 

where 

A(co) = 

In the above equation the quantities and $ii are known from the tar- 
get motion and noise characteristics as previously discussed. The transfer 
function Hf of the missile aerodynamics was chosen for this study to 
represent that for a typical variable -incidence missile. The significant 
transfer functions for this missile are given by the following conventional 
equations (ref. 7) 
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^1 %s + 1 

^ Td s(Ta^s^ + 2^g^Tg^s + l) 

^ ^ Td Vsg + 2^^T-bS + 1 
'l^ s(%s + l) 

from which the fixed network becomes 

yjl 1 + 2?i)Tb.s + 1 

- 6 ■ Tp" s2(Ta^s2 + 2^g^TaS + l) 


(T) 

( 8 ) 

(9) 


Table I svunmarizes the parameters for this particular missile at a given 
flight condition. It should be pointed out for later comparison purposes 
that these aerodynamics correspond to those used in the previous noise 
study, reference 3. The remaining parameter p is not known so that 
equation (6) can only be used to give the general form of Hco* To 
illustrate this, the form of Hco is derived in Appendix A where it is 
shown that 


2 sS (Tg^sg+2 CgToS+l) (Ta^s2+2CaTaS+D 

Hco(s) Ts (T^2g2+2^yTyS+l) (T^2g2^.2^^T^s+l) (T^2 s2^2CvT^s+1) 


Many of the factors involved in this equation are /fvinctions of the unknown 
parameter p. The exceptions are (l) the aerodynamic factors, that is, 
the gain Ts^ and the second quadratic in the numerator, and (2) the terms 
due to target motion and noise characteristics, that is, the first two 
denominator factors . The remainder of the terms can only be evaluated 
after the proper value of p has been established. This value can be . 
determined from the restriction placed on the mean-square control motion. 
The analog of equation (5) in the present case is 


poo 

6^ = . / Hco(i<*')Sco(i^) 
'-I- 


( 11 ) 


In theory, equations (lO) and (ll) can be solved for the compensating 
network, H(,q, corresponding to an arbitrary mean-square control motion, 

. In this application, however, because of the complexity of the fixed 
network and the input quantities, the computation is unwieldy. Alter- 
natively the following procedure will be used here: First calculate the 

function Hco from equation (lO) for a selected value of p. After the 
result is inserted in equation (ll) , evaluation of the integral gives a 
value of 6 . By repetition of this process for a few judiciously 
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selected values of p it is possible to determine the value of p , 
corresponding to the desired mean-square restriction £p. 

As mentioned previously, Hco does not exist, as such, in the heam- 
rider guidance system. The relationship between the filter problem being 
solved and the actual beam-rider system is illustrated by a comparison of 
figures 3 and k. From figure 3 it is apparent that if the entire guidance- 
system transfer function is represented by Yq, then 


^o “ ®co% (12) 

Thus the solution for the compensating network Hco determines the desired 
over-all transfer function Yq for the actual form of beam-rider system 
shown in figure 4 . 

Unfortunately, the theory does not give the performance of this opti- 
mum system directly in terms of the resultant error. The minimum error can 
be found, however, from the following relationship (ref. 7) 

= eip^ + 


p ^ 2 2 

= / |l - Yo(iw)| $T(tj)dco + / |Yo(ito)’| Ojgdu (I3) 


Effects of control-deflection restriction. - The procedure discussed 
in the above paragraphs has been carried out for a range of values of the 
Lagrangian multiplier p. The results of this study are presented in 
figure 5 where the errors due to target motion and noise, the total error, 
and the control deflection are given as functions of p. These curves 
illustrate the nature of the restriction of the control deflection. Large 
values of p here correspond to a large restriction and therefore small 
control deflections. It is apparent from the expression ( 3 ) that as p 
decreases the restriction on control deflection is reduced until at p = 0 
the Wiener case is reached. By cross -plotting these curves so as to elimi- 
nate the variable p which has no significance in itself, the curve of 
figure 6 is obtained. This curve is fundamental to the problem in that it 
illustrates the dependence of error on available surface deflection. Here 
the minimum obtainable error according to Wiener theory is indicated by 
the horizontal line off to the right. The curve is shown dotted in this 
region in order to indicate that these values of error would require unob- 
tainable control deflections. 

The importance of the curve of figure 6 is that it represents the theo- 
retical lower limit of error corresponding to any restriction on root-mean- 
square control motion. Each point on the curve would be achieved by a dif- 
ferent guidance-system transfer function. The interesting feature of this 
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curve is that as the surface deflection is reduced from large values, the 
minimum error increases relatively slowly over an appreciable range of 
obtainable control deflections. To determine the best operating point on 
the curve it is necessary to specify the degree to which limiting will be 
allowed to occur. Since in this linear case the control deflection is 
Gaussian, an rms control deflection of half the maximum allowable deflec- 
tion would then restrict it to within these limits about 95 percent of 
the time (ref. 8). For the specific missile being used in this study the 
maximum control deflection is about 15 which would put the desired opera- 
ting point at 7*5 as indicated in figure 6. It can be seen f rom ^f igure 6 
that for this operating point the effect of this restriction is an increase 
in error of about six feet above that of the Wiener theory. It is a some^ 
what surprising result that the increase in error is so small for such a 
great reduction in surface deflection. The significance of this is, then, 
that the optimimi Wiener result can be approached closely even- with the 
operation confined to a realistic linear range. 

Characteristics of the optimum transfer functions .- It is of interest 
to consider in greater detail the transfer functions corresponding to the 
desired operating point in figure 6. This operating point corresponds 
very closely to a value of p = lO'^ (see fig. 5) which will therefore be 
used as an example. In Appendix A the transfer functions of the compen- 
sating network, closed-loop system, and open-loop system have been derived. 
The results are given below; the numerical values of the parameters are 
summarized in table II. 

s^(Ta^sg+2^aTaS+^) 

Ty^S^+2 ^yTyS+1 

m 

T/s2+2^vT^s+1 

(15) 

Ty^S^+2^yTyS+l 

(16) 

» 

Due to Wiener Due to control Due to 

theory restriction aerodynamics 


Ta^s2+2^aTaS+l 


Hco - 


(TpS+1) (Ty2s2+2^^TyS+l) T^^B^+2^^T^S+1 


Yo = 


Ta^s2+2^^TaS+l 


° (T 3 S+I) (Ty2s2+2^yTyS+l) Tu2s2+2^^T^S+1 




s^(T^s+l) 


Tx2s2+2^xTxS+1 
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It is illuminating to trace to their origin the various factors com- 
prising these optimum transfer functions. Analysis shows that terms due 
to the basic Wiener theory also appear in the modified theory. In addi- 
tion, other terms due to the aerodynamics and the control-deflection 
restriction also appear in the solution. This situation is illustrated 
by the grouping of terms in equations (l4), (15) ^ and (l6). It can be 
said that the terms which arise from basic Wiener theory correspond in 
form to the first group of terms. However, certain constants are some- 
what altered, as can be seen from table II by a comparison of the param- 
eters given for the two theories for identical target motion and noise 
characteristics. The second group of terms in the equations indicates 
that the optimum form of the filtering term due to control-motion restric- 
tion is quadratic.^ The last group pf terms in the transfer functions can 
be traced to the aerodynamics. In fact, the numerators are identical with 
terms in the aerodynamic transfer function. 

The transfer functions are ' illustrated in figures 7(a) ^ 7(1)) ^ and 
7(c); these curves represent the optimum transfer functions in the pres- 
ence of the control-motion restriction. For comparison purposes the cor- 
responding transfer functions obtained from the Wiener theory are also 
plotted on these figures. It can be seen that, in general, the transfer 
functions for these two cases are similar at the very low frequencies and 
separate to a greater extent,^, the frequency is increased. The difference 
is primarily due to filtering terms which are required to satisfy the 
restriction on control motion. 

The added complexity has certain effects on the control motion and 
error. These effects can best be illustrated by means of the correspond- 
ing power spectra, that is, by the integrands in equations (ll) and ( 13 )- 
These spectra have been plotted in figures 8 and 9» and for comparison the 
corresponding curves for the Wiener theory are shown. It can be seen from 
figure 8 that the spectrum of the control deflection for the more complex 
system attenuates far more rapidly than does that for the Wiener theory. 

It is as a result of this fact that the modified system does not limit. 

As for the error, it is seen from figure 9 that over the important fre- 
quency range the error due to noise is not altered too seriously so that 
most of the increase in the total error is due to an increase in the tar- 
get motion component. 

Effect of noise magnitude on performance. - The preceding results were 
determined on the basis of a design noise magnitude which lay between the 
expected extremes. The noise magnitude may vary because of targets of 
different size, aspect, 'and turning rate. For this reason the error per- 
formance for the optimum transfer function has been determined and is 

^It is interesting to note that from calculations not presented here 
it has been observed that the damping ratios of these quadratics always 
lie between 0.7 and 0.8 over the entire range of control motion presented 
in figure 6. 
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shown in figure 10 (curve A) as a function of the zero frequency noise 
magnitude, N. The error performance shown is optimum only at the design 
value of noise. Shown for comparison is. curve B, which represents the 
error performance derived from the Wiener theory, corresponding to the 
same design condition (see ref, 3)* difference Between these two 

curves, then, is the increase in the error as a result of the rms restric- 
tion on the control deflection. 

As a further comparison two additional curves obtained from refer- 
ence 3 are shown. First, curve C illustrates the perfomance of a system 
(with certain saturable elements) which was optimized for the fastest pos- 
sible response to a step input in the absence of noise. The performance 
can be seen to be significantly poorer than the optimimi performance given 
as curve A. Second, the best result that could be obtained in reference 3 
is shown as cuive D in the figure. The system corresponding to this curve 
was not optimum because the linear Wiener theory used in this approach was 
not capable of evaluating all of the various filtering terms discussed in 
preceding paragraphs. For example, terms due to Wiener theory were not 
altered, and terms due to aerodynamics were not used. Terms due to con- 
trol restriction were included but the performance was limited by the 
arbitrarily assumed form. Thus the system was not optimum and limiting 
occurred; the performance deteriorated accordingly. 

System synthesis .- To achieve the results which have been presented 
it is necessary to design the guidance system shown in figure 4 to have 
the optimimi transfer function Yq. This can be done in many different 
ways since the actual beam-rider system consists of two distinct parts: 
the tracking radar and the missile-control system. For example, the opti- 
mum transfer function Yq might be split up and apportioned between these 
two parts of the guidance system. Or, on the other hand, the missile- 
control system alone might be designed to approximate the optimum Yq in 
which case the tracking radar should be designed to have a relatively fast 
response. In a similar manner the tracking radar could be optimized, in 
which case the fast response should be designed into the missile-control 
system. 

To investigate ‘all these possibilities is beyond the intended scope 
of this report. As an example, however, of one possible design the choice 
of missile -control- system filtering was investigated. For this case it is 
desired to design the missile-control system to match the optimum transfer 
function Yq. This procedure is not unique; the analysis used herein is 
based, primarily on cut-and-try procedures. It has been found that the 
desired system could be synthesized in a conventional fashion illustrated 
in figure 11. To show how this system is capable of supplying the optimum 
filtering the open-loop transfer function has been derived in Appendix B 
where it is shown that with certain assixmptions 
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Yu _ kgkgV (Tips + 1)(T4S + l)(Tb^s^ -I- 2^^)TbS + l) 

^ I Td + ks ^3 s^(ThS + l) (TgS + l)(as® + + cs + l) 

(17) 

Comparison with the optimum Hp in equation (l6) shows that certain dif- 
ferences in form exist. However, equation (l6) can he expressed in another 
form hy factoring approximately certain quadratic terms which have large 
damping ratios. Equation (l6) can he written then as 

(TgS + l)(TgS + l)(Tb^s^ -t- + l) 

^ b 2(T^S + l)(TxS + l)[(TxS + l)(Ty2s2 + 2^yTyS + 1)] 

( 18 ) 

It is now possible to choose the parameters in equation (I 7 ) to match those 
of equation (I 8 ); the specific values are given in table III. 

The remainder of the system consists of the tracking radar whose 
response should he made relatively fast for this design approach. Further, 
from the plots presented in figures 8 and 9 H is apparent that the track- 
ing radar should he designed so as not to alter appreciably the spectral 
distribution of error or control motions. Thus a break point or natural 
frequency of several cycles per second would be satisfactory. 

Comparisons and other considerations.- In a REAC simulation of this 
system its operation was examined in further detail. Of particular impor- 
tance in this examination was the control deflection and the associated 
servo energy. It was found that the control deflections were effectively 
restrained to within the linear range and as a result the servo energy was 
greatly reduced. For a given servo the power expended is proportional to 
the time average of the sum of the absolute displacements of control 
motion. Thus the average servo power over a sufficient time interval of 
duration t can be determined from 


Lj OU I 

average servo power — - — 

For this system there is more than a two-thirds reduction in required servo 
energy over that for the system with smallest servo energy requirement dis- 
cussed in reference 3, The saving in servo energy becomes even greater 
when compared to systems for which the error performance becomes progres- 
sively worse than the optimum. The control-deflection time history 
obtained from the REAC was used to show that the required control rates 
are easily attainable since they rarely exceeded 3 "to 3*5 radians/sec. 

This system has other virtues. For example, there is .also a reduction 
of voltages within the circuit to reasonable and easily obtainable values. 
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From this simulation work, then, two facts are apparent. First, it 
is clear that the validity of applying only a control-motion restriction 
in the analysis has heen verified, since the control rate and system vol- 
tages did not saturate. Second, it is apparent that in comparison with 
systems of reference 3 which limiting occurred, the present system not 
only produced a smaller error hut was accompanied by other desirable 
effects such as decreased servo energy and circuit voltage requirements. 

The missile-control system discussed in the above paragraphs is 
intended only as an example designed to achieve optimum performance 
against noise. It is clear that since the system operation is essen- 
tially linear, all designs which approximate the desired . Yq would give 
approximately the same performance. However, it is often necessary to 
satisfy other requirements which are related to the specific design objec- 
tives of' the system and which may favor certain designs or system con- 
figurations. Because of the freedom in choice in distributing the .opti- 
mum filtering in the system, the optimijm characteristics specified by the 
theory do not, in general, place an inherent limitation on satisfying other 
requirements. For example, requirements on transient response during the 
beam entry phase may dictate a design in which more of the optimum filter- 
ing is located in the tracking radar. Similarly, requirements of simpli- 
city on certain parts of the system may dictate yet a different design. 
Because these and similar requirements are related to the detailed design 
objectives, an investigation of such factors is beyond the scope of this 
report. 


CONCLUDING REMAEKS 


This study has considered the application of Newton^ s modification 
of the Wiener filter theory to the choice of optimum transfer function 
for the beam-rider guidance system^. The analysis has been applied to a 
typical variable -incidence missile at a given flight condition attacking 
a maneuvering target in the presence of glint noise. By minimizing the 
mean-square error with a side restriction on the mean-square value of 
certain quantities, limiting effects, of primary importance in missile 
guidance, are largely eliminated. 

Although many limiting-type nonlinearities are of importance, this 
study has shown that the critical nonlinearity is due to control-deflection 
limiting, since a restriction applied to the Reflection sufficient to con- 
fine control motions to within the linear range also satisfactorily pre- 
vented other types of limiting. Results have been given which illustrate 
the effect of the control -motion restriction on the minimum rms error. 

Zero restriction corresponds to the Wiener theory whose results are 
unattainable since impossibly large control motions are demanded. How- 
ever, as the surface deflection is reduced by means of an increasing 
restriction, the minimum error increases slowly until at a realistic 
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value of the control deflection the error has increased only a few feet. 
Thus it is clear that the optimum performance predicted hy the Wiener 
theory can still he approached closely hy confining operation of the sys- 
tern to a realistic linear range. 

The transfer functions of the optimum guidance system have heen given 
and have heen compared to the corresponding transfer functions from the 
Wiener theory. Differences were shown to he the result of added filtering 
terms demanded hy the restriction on control motion. The various filtering 
factors can he traced to their origin. It was shown that terms due to 
basic Wiener theory appear in the solution, and that the additional terms 
can he associated with the control-motion restriction and the missile aero- 
dynamics . 

For the optimum system, performance against variable noise magnitude 
was given and was compared to Wiener theory. As an example of one possible 
design 'which achieves these results^ a missile-control system was synthe- 
sized to match the optimum transfer function. In a REAC simulation the 
operation of this sytem was examined in further detail. It was found that 
the magnitude of circuit voltages, and control-motion rates are greatly 
reduced in comparison with systems in which limiting occurs. Similarly, 
the required servo energies are appreciably reduced, resulting in about 
a two -thirds saving. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 11, 1955 
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APPEMDIX A 


DETERMINATION OF THE OPTIMUM MODIFIED WIENER FILTER 


The solution for the optimum compensating network Hqq illustrated 
in figure 3 is given "by equation (6) in the text as follows: 




-iu)t r“ Hf(io-)<I>T(oc)e 
A" (a) 


iat 


da dt 


where 


A(co) = 


Hf (io))Hf (iw) + p <I>ii(w) 




> (Al) 


In this equation $11 is defined as 


<l>ii = 


(A2) 


where and <Iijj represent the spectral density of the target motion and 
noise, respectively. The quantities A'*’ and A~ are defined as the factors 
of A with poles and zeros in the upper and lower half -planes, respec- 
tively. Thus 


A = A'^A" 


(A3) 


For the case to he considered here the target motion and noise will he 
taken as 


K0l)^((j^ + k^) 

% = N (A5) 

It will he noted that the noise spectrum has been approximated here by a 
constant in order to reduce the complexity of the calculations. As shown 
in reference 3, this approximation is valid because the noise defined by 
equation (l) is essentially flat compared to the passbarid of the optimum 
system. It should also be pointed out that use of equation (A4) leads to 
certain mathematical difficulties which occur in the process of evaluating 
the right-hand side of equation (Al) , because the theory requires that the 
poles of not be located on the real axis. To avoid these difficulties 

0 
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it is necessary to modify the target-motion spectral density to the 
following 

+ k2) 

vhere t)^ arid are any small real numbers. The solution to equation 
(A1) vill then he a function of iw, and t)^. The desired answer is 
obtained hy taking lim Yq as and t|^ 0. However, it can 

he shown (see ref. 3) that the same answer can he obtained more simply hy 
taking and equal to zero as would he obtained from the rigorous 
process described above. 

Now hy combining equations (A2) , (aU) , and (A5 ) , 


^ii(°') = 

- M + (ka^/ri?j) 

a'*(a^ + k^) 


(g - api)(g - apg) (a ~ apa) (a + api) (a + apa) + -^ps) 
g^(g + ik) (g - ik) 


g^(g + ik) (g - ik) 

where the gpm represent ■ roots in the upper half -plane. 
p(a) is 

3 

p(g) = n (g - gpm) 
m=i ^ 

s g® + hga^ + "big + ho 


(AT) 

The polynomial 


(a8) 


where 


h2 “ “(^pi + ^p2 ^ps) 

hi = ^('PlCtpg + o^piOt-ps + ^p2^p3 

ho = “Qt'pio^pa^ps 


(A9) 
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In a similar manner 


-p(-a) = a® - bgct® + bia - bo 


(AlO) 


An expression to be used later can be obtained from these preceding 
equations as follows 

a® + k®a'* + = (a® + ba^® + bitt + bo) (cx® - b 2 a® + bia bo) 


from which 


2 _ 


(All) 


The quantity Hf in equation (A1) represents the fixed network or 
aerodynamics and its general form (ref. 7) is given by 


nWia) -M- 1 t(i - . l(2g„T^»)l 

Hfua) - 6 - - ^ KSSaTaa)] 


(A12) 


and 


[(1 - T^,®g®) - i(2g^T^g)] 

Hf(ia) - - ^^2 a^[(i - Ta^a^) - i(2LT„a)] 


(A13) 


Thus 


5-T^tt N - 1 [Tb^g^ + (^^b%^ - 2Tb®)g® + 1] 

Hf(ia)Hf(ia) a^[Ta^a4 + (4^^%^ “ 2Ta®)a® + 1] 


(All^) 


Accordingly it can be shown that in equation (Al) 


Hf(ia)Hf(ia) + p 

^ pTs*Ta*g^->-pTs^(4^a^Ta®-2Ta®)a^+(Tb^+pTs^)a^+(Ugb^Tb®-2Tb®)a®+l 
Tg-^w^ [ Ta‘*g‘*+( 4?a^Ta® -2Ta® ) a®+l ] 


n (g-gqn) n (<x+aqn) 


■ a*tTa‘a*.(!l?a=Ta=-2Ta=)a=.ll 


n=i 


= pTa^ — 


q(a)q(-a) 


a-* [ Ta-^a^+(4^a%^-2Ta^) 


(A15) 
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where the aq^ represent roots in the upper half -plane. 
q(a) above is 

q(a) = n (a - aqn) 


The polynomial 


n=i 


= + cqpP + + Citt + Co 


(ai6) 


where 


cs = - (o.qi + «qa + ^qs + aq4) 

Cg = OLqiOCqg + Q^qiO^qs + ^^i^q4 ^qs^Q3 ^qs^(j4 ^qs^^4 

Cl = ” (ctqictqaaqs + °<3i°<j2Cq4 + °^qi*^q3Ctq4 + ctq2Ctq30q4) 

Cq = ciqia.q2aq3a<}4 
In a similar manner, 

q(-a) = a'^ - caa® + CgO.^ - Cia + cq 


> (Al?) 


(A18) 


An expression to he used later can he obtained from these preceding equa- 
tions as follows 

\ 

pTg^ Ta'* aS+pTg^ -2Ta^ ) o-®+ ( V +pTs* ) ^h^ V -2T^® ) a^+l 


= pTg^Ta^q(a)q(-a) 

\ 

= (fPs'^Ta^ (o^+cso®+C2a®+Cia+Co) (c^ -C3a®+C2a®-Cia+Co) 


from which 



(A19) 


The second equation in (Al) can now he found by combining equations (A7) 
and (AI 5 ) 

^ q(a)q(-a)p(a)[-p(-g)] 

A ~ a®[Ta‘*g'* + “ 2Tg^®)a^ + 1] (a + ik)(a - ik) 

(A20) 
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This expression can be split into two factors with poles and zeros in the 
upper and lower half -planes. From the preceding definitions q(a) and 
p(a) both have their zeros in the upper half -plane. , Thus 


A+(a) 


pTa“^N 


q(a)p(.g) 

a-*(l - Ta^a^ + i2^a,'Tad) (a - ik) 


(A21) 


q(-a)[-p(-g)] 

a^(l - Tg^^a^ - i2^a.'^a.a) (a + ik) 


(A22) 


Now, let us consider the evaluation of the first integral in equa- 
tion (A1) . This integral can be evaluated by combining equations (A4) , 
(AI 3 ) , and (A22) as follows: 




1 


r°° Hf {a)$T(a)e^®'^ 


A~(g) 


da 


ka^ 

Ts^it 



(1 - T^Sg^ - i2g^T^g)e^”'^ 
g2(a - ik)q(-a)[-p(-g)] 


da 


ka^ 



(1 - - i2^T3T^g)e^^^ 

4 3 

a^(a - ik) IT (a -h aqn) H (a + apm) 
n=i iu=i 


da 


2 

= - I f(g)<ia (A23) 

Tg2« J_^ 

where f(a) is defined as the integrand. By considering a to be a com- 
plex variable, the integral can be evaluated by a contour integration indi 
cated by the following sketch: 
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There are two poles within this contour, a second-order pole at the' origin 
and a simple pole at a = ik, as indicated. The residues are 


Res(o) = lim [a%(a)] 
a-»o da ' 


t khoCo2^-bTb + ikCbpCx + biCp) + hpCo 

kho^o k^ho^^o^ 


— TtT 

khoCo 


(A2U) 


Res(ik) = lim (a - ik)f(a) 
a-4ik 

^ (Tb^k^ + 2g^T|,k + l)e~^^ 

-k^q(-ik)[-p(-ik)] 

= (A25) 

Now, hy means of the residue theory, the value of the integral is given 

■by 


[Res(o) + Res(ik)] 
•*-s 


Thus from equations (A24) and (A25) ^ 




2ika^ 

Ts^ 



+ 


t 

khoCo 


+ 72® 



(A26) 


The second integral in equation (Al) is denoted by I 2 and is merely 
the Fourier transform of Ii. Hence 



'1 kboCo '2 J 

(A27) 


2ika^ r“ 

Ts^ 1 
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Thus it can be shown that 


2ika^ (7i + 72)'boCo(ia))^ + [kritpCp + (l/k)]la) + 1 
Ts^'boCo (iw)®(k + iw) 


21ka^ + 2^a.Ta(ito) + 1 

Ts^boCo (io))^(k + iw) 


(A28) 


It should he noted that since the r\^s were assumed to he zero, questions 
concerning the existence of the above integral arise. However, as indi- 
cated in reference 3^ T^en the analysis is made without this assumption, 
the integral in equation (A2T) does exists and the limit of this integral 
as the Ti*s approach zero becomes precisely equation (A 28 ) . 

The coefficients can he expressed in a more convenient form by elimi- 
nating the intermediate parameters and y^ by means of definitions 
given in equations (A2ii-) and (A25) . The following results are then 
obtained: 

= k7i^oCo + ^ 




^ 0^1 + ^ 1^0 
boQo 


(A29) 


Also 

Ta^ = (7i + 


_ , kl3oCo2^-^T-^j + ik(boCi + ~biCo) 1 

k^q(-lk) [ -p(-ik) ] k^bpCo k^ 

(A30) 

For the missile parameters, of table I and the value of k given on page 28 

. + 2^-[jTi3k + 1 « 1 


so that very nearly 


Ta^ = 


^qCq 


. L 

k^q( -ik) [ -p(-ik) ] k k^ 


(A31) 
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Now, from equation (Al) the optimum compensating network can he 
determined as follows: 


Hco(iw) = 


I2 


(A32) 


2jtA‘*’(w) 

Substituting equations (A2l) and (A28) into (A32) gives 
Hco(i*^) 

( io)) ^ ( io)) ^+2 iai) +1 ] [ ( iu) ^+2 (i(j^) +1 ] 

rf^pTs^Ta^hoCo p((j)q(u)) 

(A33) 

To reduce this expression to a more standard form it can he seen from 
equation (A16) that 

q(w) = “ “qn) 


= Co n (—•+ 1 

^ n=i V ^'•qn 


Similarly from equation (a 8) 

3 

P(‘^) ^ mSi 

= ho n (1 ^ + i") 

m=i \ a.pm J 


(A3»^) 


- (A35) 


Also hy using the definitions given in equations (All) and (AI9) the com- 
pensating network can he shown to reduce to 

2co(iu)) = Ts^ (i(^)^[Ta^^(ia))^ + 2^g_To^(iuj) + 1 ] [ ( io)) ^ + 2ga,Ta(i^) + 1] 

^ ^ + 1) 

m=i V / n=i \ J 

(A36) 

By the substitution of the conventional complex frequency s for iio 


Hco(s) - Tg 


s^(T^^s^ + 2^„T„s + 1) (T,,^s^ + 2^„T,,s + l) 


(A3T) 
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It is shown later in the numerical examples that certain terms combine to 
give quadratic factors as shown helow 


Hco(s) - Ts^ 


where 


s^ (T/s^+2gQ^Tg,s+l) (Ta^s^4-2^a'TaS-»-l) 


(Tps+l) (T/s^+a^yT^s+l) (Tu 2 s 24 fi^^T^s+l) (Tv^sS+a^^TyS+l) 

(A38) 


To - ^ 

rp ■ 2 

1 

^ P api 

. Tu = - 

oqiOqa 

= . 1 
7 Ctp2dp3 

2 ?uTu = 1 

oqi + aqa 
otqiaqa 

i ^P 2 + O'PS 
^7 7 ttpattpa 

1 

II 

1 

oq3aq4 


2 ?vTv = i 

aqs + aq4 
ocq3aq4 

The compensating network given 
the over-all transfer function 

in equation (A 38 ) 
Yq. From figure 

can be used 
2 


Yo(s) = Hco(s)Hf(s) (A39) 

Substitution of equations (A12) and (A 38 ) into (A39) gives 


Yo(s) = 


(T„.^s^+2^^T^s+l) (T-h^s^+2^^T^s+l) ^ 


(Tps+l) (Ty2s2+2?yTyS+l) (T^2 s2^2^^T^s+1) (T/s2+2^^T^s+1) 


(AlK)) 


If Yo is considered to represent a unity feedback system, the equivalent 
open-loop transfer function is shown in the later specific example to be 
of the following form: 


^lo(s) = 


Yq(s) 

1 - Yo(s) 


^ + 2CaTgS + l)(T-h^s^ + 2^^T^S.-h 1) 

^ s^Ct^^s + 1)(Tx^s 2 + 2^x^x®+in'^y^s^ + ^^y^y® 

Uhl) 
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Because and are close enough to unity (table II) it is possible 
to factor the corresponding terms approximately to give 

,,(S) . k, (^S 1) (T^s l) (Tb^s^ ■. 2^^TbS ^ 1) 

s2(Tp^S + l)(TxS + l)[(TxS + l)(Ty2s2 + 2^yTyS + l) ] 

(Ai<-2) 

The numerical evaluation of the optimum system has been carried out 
for the following values of target motion and noise: 

k = 0.4 a = 1 g 

T = 5 sec W = 15 ft^/radian/sec 

From this example p will be taken to be 10“* since, as shown in the text 
this value results in a near-optimum restriction on control deflection. 

Evaluation of the numerator in equation (AT) gives 

a® + k^a'* + = a® + O.16 a'^ + 8.8 

3 3 

m=i ‘ mSi 


where • 

api = i 1.456 

ttp2 = 1.427 exp(i 0.543) 

ttp3 = -1.427 exp(-i 0.543) 

which are all located in the upper half of the complex a plane. Then 
from equation (A9) 

hs = -i 2.94 
bi = -4.201 
bo = i 2.97 
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Similarly from equation (AI5) for p = 10'*, 

pTs%*a®+pTs*(4^a^Ta^-2Ta^)a®+(Tb'*+pTs^)a*+(U?T32^2.2T^2)a2+i 

= (22. 58xlO"®)a®-(7A75XlO'S)a®+(6. 267 x 10 ■®)a*-(0.00609l^)a2+l 

4 4 

= pTg^Ta** (a-oqn) (a-HXqn) 


Solving for the roots gives 

otqi = 3*593 exp(0.804) 

Oqa = -3*593 exp(-0.80U) 

Oq3 = 12.92 exp(0.80U) 

ocq^ = -12.92 exp(-0.80l|) 

which again are all in the upper-half a plane. From the definitions in 
equation (AI7 ) } 

C3 = -i 6.677 
C2 = -187.61 
Cl = i 883*1^0 

Co = 2155 


Now from the values of the roots apm and oqn 'well as the h and 
coefficients', the optimum compensating network can he calculated by means 
of equation (A38) . For this purpose it is necessary to evaluate the left- 
hand sides of equations (A29), (A3I), and the parameters following equa- 
tion (A38) . Results of these calculations are given in table II in terms 
of damping ratios and undamped natural periods. Substitution of these 
parameters in equations (A38) and (a 40) gives 
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APPENDIX B 

MISSILE-CONTROL-SYSTEM APPROXIMATION TO 


THE OPTIMUM TRANSFER FUNCTION 


The optimum linear system is illustrated in figure 11. In the follow 
ing section the system equations are derived. 

Using figure 11 it can. he shown that 

^ ^ ks(TMS + l)(TaS + l) 

% T^s(Tj^s + l)(Ta^s2 + 2^a'^aS + 1)(T3S + l) + kekgsCTj^s + l) (TgS + l) 

(Bl) 

To simplify equation (Bl) it is convenient to make T 3 = T(^. This assump- 
tion is not essential hut its use leads to simpler equations. In certain 
cases where gust disturbances are serious it may he more desirable to 
choose T 3 small. With the former choice 

_ kg(TMS+l) 

En " s[(TeLTa2T^)s3+(T^Ta2+T,i2^aTaTi)s^+(TdTi+T<i2?aTa+ksk3T2)s+(T^+ksk3)] 


ks (%s+l) 

Td+kgk 3 s(as®+hs^+cs+l) 


where 


Td + kg kg 


h = 


Td(T^^ + 2^^T^T,) 


Td + kgkg 


^ ^^a'^a'^d '^I'^d ^s^a'^a 

Td + ksk3 


(B2) 


(B3) 
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Thus the entire open-loop transfer function can he written as 

^ ^ ^ kgkgV (TiqS + 1)(T4S + 1)(T-q^s^ + 2^^TbS ■+ l) 

T(i + kgkg b^(TiiS + 1)(T5S + l) (as® + hs^ + cs + l) 

(b4) 


By comparison with equation (Ah2) it can he seen to he of the same form as 
the optimum transfer function. Thus the following correspondences between 
quantities in equations (B^J^) and (Ah2) are appropriate: 

Til ^ T>, 

Ts ^ Tx 
Tio Ta 
T4 Ta 

•/ 

Also it is apparent that 


a » TxTy^ 

h -e^TyS + Tx2^yTy > 
C «-» Tx + 2^yTy 


(b6) 


For these values of a, h, and c, equations (B3) can he solved for the 
remaining parameters to give • ^ 


Ti 


aTa 

- a2^a 


kgks 


T(jTa®Ti - aTd 
a 


> (B7) 


T2 


cT^j^ + ckgkg - g^g^Tg^T^j^ - Tj^Tj^ 
ksks 
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The value of kgks is determined from the desired gain: 

Tf5 + kek« ' , 

k2ks <— > k^ — — — (b8) 

Since only kgkg and kska are specified, one of the three gains kg, ks, 
or ks may "be chosen arbitrarily, subject to the condition that voltage 
limiting does not occur. From the above equations the system parameters 
have been determined and are tabulated in table III. 
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TABLE I.- SUMMARY OF AERODYNAMIC PARAMETERS FOR EXAMPLE MISSILE 


Parameter 

Value 

Ta 

0.0775 

Tb 

.0552 

Td 

2.087 

Ts® 

.0007911 

% 

. 81^6 


.0536 


.0220 


TABLE II . - SUMMARY OF PARAMETERS OF OPTIMUM TRANSFER FUNCTIONS 


Parameter 

Modified Wiener 
tlieory value 

Wiener theory 
value 

Ta 

1.13 

0.925 


.805 

.765 

Tp 

.687 

.687 

Ty 

.700 

.700 


.519 

.519 


2.5 

2.5 


3.2 

7.k2 

Tu 

.278 

— 


.720' 

— 

. Tv 

.0774 

— 


.0581 

— 

Tx ‘ 

.184 

— 


.722 

— 

Ty 

.0787 

— 

^y 

.0513 

— 
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TABLE III.- SUMMARY OF PARAMETERS OF OPTIMUM SYSTEM 


Parameter 

Value 


0.0538 

T4 

1.13 

Tio 

1.13 

Til 

2.5 

T5 

.184 

ks 

.0463 

Ti 

.186 

ka 

-.674 

Ta 

.373 

T3 

.846 
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Radar beam 



Figure 1.- Beam-rider guidance system. 







Fixed network 



Figure 2*- Block diagram of general filter problem with restrictions 
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Figure 3 .- Block diagram of filter problem for the beam-rider system. 


Missile-control system 



Figure 4.- The "beam-rider system. 
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(a) Compensating network . 

CO 

Figure 7*" Optimum transfer functions. 
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Frequency.eps 

(b) Closed -loop transfer function Y^. 

Figure 7*~ Continued. 
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Figure 11.- System approximation of optimum transfer function with control motion restriction. 








